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We present first results of a study of two flavour QCD with Wilson fermions at flnite temperature. We have 
used tree level Symanzik improvement in both the gauge and fermion part of the action. In a first step we explore 
the phase diagramm on an 8'^ x 4 lattice, with particular emphasis on checking Aoki's conjecture with an improved 
action. 



1. Introduction 

The study of the finite temperature phase tran- 
sition in QCD with Wilson fermions is much more 
complicated than in the staggered fermion for- 
mulation, because of the absence of an order pa- 
rameter due to explicit chiral symmetry breaking. 
This means e.g. that the existence of a massless 
pion at finite lattice spacing is not at all obvious. 

In recent years a detailed picture of the finite 
temperature phase diagramm of QCD with Wil- 
son fermions has emerged [1]. This picture is 
based on the idea of spontaneous breakdown of 
parity and flavour symmetry [2] and has been 
investigated analytically as well as numerically. 
The features of this phase diagramm are (i) the 
critical line KciP) defined by a vanishing pion 
screening mass for finite temporal lattice size 
turns back towards strong coupling forming a 
cusp, (ii) the region bounded by the critical line 
represent a phase of spontaneously broken parity 
and flavour symmetry, (iii) the finite temperature 
phase transition line K,t{P) presumably does not 
cross the critical line, but runs past it towards 
larger values of the hopping parameter [3]. 

From an analysis of the Gross-Neveu model in 
two dimensions, where three cusps connected to 
doublers develop, one expects the critical line for 
QCD in four dimensions to form five cusps mov- 
ing towards weak coupling with increasing tempo- 
ral lattice size Nr ■ Simulations with the standard 
Wilson formulation for quarks and gluons have 



*Work supported by NATO grant no. CRG940451 



shown, that at N^- = \/{aT) = 4 the tip of the 
cusp lies in the strong coupling regime at « 3.9 
and moves only slowly towards weak coupling as 
Nr is increased. 

Since one expects the same features to hold for 
a wider class of actions including the clover ac- 
tion, which tends to reduce cutoff dependencies, 
a study of the phase diagramm using improved 
actions is of practical as well as theoretical im- 
portance [4]. 

2. Results 

We have conducted a simulation of 2 flavour 
QCD on an 8'^ X 4 lattice using tree level Symanzik 
improved actions for quarks and gluons. In the 
gauge sector this amounts to adding a 2 x 1- 
loop to the standard plaquette action and for the 
fermions in adding the clover term with csw = 
1. We have used a Hybrid Monte Carlo algo- 
rithm with a timestep 5t = 0.01 and the num- 
ber of molecular dynamics steps Nmd — 20, 
which so far amounts to rather short trajecto- 
ries. We simulated several k. values for each 
(3 = 3.00, 3.50, 3.75 and 4.00. We have measured 
the Polyakov loop, the pion norm and the average 
number of iterations it takes to invert the fermion 
matrix. Each observable is now discussed in de- 
tail. 

2.1. Locating the critical line 

Although not a physical observable in its own 
right, the average number of iterations it takes to 
invert the fermion matrix is a very good indica- 
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Fig. 1 Average number of iterations as a function of k 
for P = 3.00 (triangles), /3 = 3.50 (circles), (3 = 3.75 
(diamonds) and = 4.00 (boxes). The inverter used 
are conjugate gradient (lines), overrelaxed minimal 
residual (dashed), BiCGstabl (dotted) 



tor for criticality. The simulations where done in 
the following way. At each /3-value we started a 
simulation at k = 0.12 and used a thcrmalizcd 
configuration from this run as a start configura- 
tion at a higher value of k. We continued to do so 
for higher and higher K-values. At those /3- values, 
where we saw a drastic increase in the number 
of iterations, we started a simulation at a much 
higher K-valuc and continued towards smaller re- 
values. Our findings are summarized in Fig. 1. 
At /? = 3.00 we were not able to simulate the 
system for 0.1770 < k < 0.1825. At /3 = 3.50 
we saw an initial increase of iterations, even af- 
ter switching from minimal residual to conjugate 
gradient, which usually decreased the number of 
iterations. With increasing k the number of it- 
erations decreased again, only to increase again 
at fairly high k- values. At /? = 3.75 and 4.00 
we saw a similar behaviour as at /3 = 3.50 but 
not as pronounced. We experimented with dif- 
ferent inversion routines and our conclusions are, 
that close to the critical line conjugate gradient 
is superiour to overrelaxed minimal residual and 
BiCGstabl. 



2.2. Pion Norm 

Since it is not possible to reliably extract the 
pion screening mass on small lattices, wc use the 
pion norm instead to indicate the existence of a 
critical line of vanishing pion screening mass. The 
pion norm is the integrated pion correlator and is 
defined as follows: 
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Here A4 is the fermion matrix on a particular 
gauge configuration. Near the critical line the 
pion norm behaves as 11 « 1/n^J, hence a di- 
verging pion norm indicates the existence of the 
critical line. Our results are displayed in Fig. 2. 
At /3 = 3.00 we find a clear signal for two critical 
lines close to n = 0.1770 and 0.1825. The differ- 
ence in K is already quite small, so we are near the 
tip of the cusp. At /3 = 3.50 the pion norm devel- 
ops a small peak at k = 0.1550, which is located 
where the crossover from the low temperature to 
the high temperature phase starts (see section on 
Polyakov loop). No divergent behaviour is ob- 
served and the critical line ceases to exist in this 
coupling region. The simulations at /3 = 3.75 and 
4.00 do not even find a peak for the pion norm, 
but rather a smooth behaviour as a function of 
hopping parameter. We conclude, that the crit- 
ical line turns back towards strong coupling and 
develops a cusp between /3 = 3.00 and 3.50. 

2.3. Polyakov Loop 

The Polyakov loop is defined as follows: 

^=]^E]^Tr f[U,ix,r) (2) 

This observable is sensitive to the finite temper- 
ature phase transition although it is no order pa- 
rameter in the full theory. Our results are dis- 
played in Fig. 3. At /3 = 3.00 we find a confined 
phase for k < 0.1770. For for k > 0.1825, when 
one approaches the critical line from above, the 
Polyakov loop decreases to \L\ « 0.1. This indi- 
cates that the system develops a low temperature 
behaviour in the vicinity of the critical line. On 
the other hand we do not see a sharp crossover, so 
we cannot conclude that one crosses the thermal 
line as one lowers k towards Kc{(3). 
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Fig. 2 Pion norm as a function of /t for /3 = 3.00 (tri- 
angles), p = 3.50 (circles), /3 = 3.75 (diamonds) and 
/3 = 4.00 (boxes). 



Fig. 3 Polyakov loop as a function of k for /? = 3.00 
(triangles), /? = 3.50 (circles), (3 = 3.75 (diamonds) 
and (5 = 4.00 (boxes). 



At /? = 3.50 the Polyakov loop displays a sharp 
crossover phenomenon, which means that the sys- 
tem crosses the thermal Hne for 0.1550 < k < 
0.1600. 

At /3 = 3.75 and 4.00 the system is in the high 
temperature phase down to k = 0.12. This is not 
unexpected, since the finite temperature phase 
transition in the quenched theory for our choice 
of action occurs at /3c = 4.07. 

3. Summary and Conclusions 

From measuring the pion norm, Polyakov loop 
and the average number of iterations to invert the 
fermion matrix we conlude, that the finite tem- 
perature phase structure observed with the stan- 
dard Wilson formulation is preserved for the tree 
level Symanzik improved formulation of quarks 
and gluons. We note that the difference in f3 be- 
tween the location of the cusp and the quenched 
/?c is considerably reduced. What this means 
in terms of physical scales remains to be inves- 
tigated by measurements of the lattice spacing. 
After this preparatory simulation on a small lat- 
tice, we will investigate the phase diagramm on 
larger lattices including a precise measurement of 
the pion screening mass and quark mass, the lat- 
ter enabling us to study the chiral condensate as 



well [5]. 
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